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We settle the equivalence between the problem of hitting a polyhedral set by 
the orbit of a linear map and the intersection of a regular language and a language 
of permutations of binary words (Pi-realizability problem). 

The decidability of the both problems is presently unknown and the first one is 
a straightforward generalization of the famous Skolem problem and the nonnega- 
tivity problem in the theory of linear recurrent sequences. 

To show a 'borderline' status of JVrealizability problem with respect to com- 
putability we present some decidable and undecidable problems closely related to 
it. 

This paper is an extended version of the journal publication [16] and contains 
some additional results. 



Introduction 

Let $ be a linear map of a vector space V into itself and let x G V be a vector fin V. 
The iterations of <f> applied to x define an orbit Orb$ x, i.e. the set 

{<S> k x : k € Z+}. 

In the present paper we discuss algorithmic issues related to orbits. We assume 
that V is a rational coordinate space, so 4> and x are represented by their (rational) 
components. 

An orbit description problem consists in finding specific relations which either hold 
for all vectors in the orbit or are violated by at least one vector in the orbit. Here we 
limit ourself to a simple case where the relations are formed from Boolean combina- 
tions of linear equalities and inequalities (see the exact definition below in Section[2]). 

Note that the most important case of the orbit problem is the chamber hitting prob- 
lem. In this particular case we check whether an orbit intersects a closed polyhedron 
(i.e. a solution set for a finite system of nonstrict linear inequalities). 
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The orbit description problems are related to some problems on linear recurrent 
sequences. 

A linear recurrent sequence (LRS) x n of a degree d is defined by: 

when n > d, 

(1) 

when 1 < n < d, where a^, bj are constants. 

The famous Skolem problem is perhaps the most known algorithmic problem on 
linear recurrences. 

The Skolem problem. Let {x n } be a specified LRS with integer coefficients. 
Whether Xk = for some kl 

The decidability of the Skolem problem is presently an open question, but it is 
known that it is decidable for the degrees < 5 (the cases d = 3, 4 are worked out by 
N. Vereshchagin [ 17 1, and the case d — 5 is solved by V. Halava et al. Q). 

In the opposite direction, the best hardness result on the Skolem problem is NP- 
hardness 0. 

LRS is called nonnegative if all its elements are nonnegative. One more important 
algorithmic question about LRS is called the nonnegativity problem and consists in 
checking nonnegativity of a LRS. This problem is at least as hard as the Skolem prob- 
lem. Being a bit more formal it means that the Skolem problem is Turing reducible 
to the nonnegativity problem. This statement certainly belongs to the public mind, but 
see the proof of the statement in Section[4] 

The nonnegativity problem is decidable for LRS of the degree < 3 (V. Laohakosol, 
R Tangsupphathawat [9]). 

In this paper we often use some standard constructions from the theory of algo- 
rithms and the complexity theory, e.g., Turing reducibility, m-reducibility, polynomial 
reducibility. In particular, Turing reducibility of a Problem I to Problem II means that 
while solving Problem I the reduction may ask an oracle, who can give answers to 
Problem II (see, further Ifl3ll6l.) 

The relations between the orbit description problems and LRS are described in 
Section[2] 

These results are certainly not new but we include them for completeness sake and 
in order to introduce necessary notation and terminology. 
Let briefly sketch the contents of the paper. 

We recall basic properties of LRS in Section Q] LRS are closely related to regular 
languages. In particular, any LRS can be represented as a difference of the generating 
functions of a pair of regular languages (see, e.g., ifTTl Cor. 8.2]. We will use a modifi- 
cation of this result (see Theorem[3]from Section|4|. (All necessary facts about regular 
languages could be found in (Sj [2j .) 

The main result of the present paper consists in an algorithmic equivalence be- 
tween the chamber hitting problem and checking a particular property of the regular 
languages. Namely, the property involved consists in checking whether a regular lan- 
guage contains at least one word from a special set of words. We call this set a permu- 
tation filter. Speaking informally, an arbitrary word from the permutation filter gives a 
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permutation of all binary words of a fixed length n. See a formal statement in Section[3] 
below. 

The permutation filter is denoted by Pg and the corresponding problem of checking 
this property of regular languages is called Pn-realizability problem (see, Section[3]l. 

The proof of an algorithmic equivalence of a problem of Pg-realizability and the 
chamber hitting problem (Theorem[2]from Section[3]) proceeds in several steps. 

At first, we describe in Section [4] a polynomial reducibility of the chamber hitting 
problem to the Pn-realizability problem. The reducibility uses the aforementioned 
Theorem[3] 

As a consequence of this result we prove (see, theorem [4]i NP-hardness of a Pb- 
realizability problem. This result may be of independent interest. 

To construct a reduction in the opposite direction, i.e. a reduction of the Pr- 
realizability problem to the chamber hitting problem, we settle some technical difficul- 
ties. It turns out that a natural construction described in Subsection |5. II gives a reduc- 
tion of the PB-realizability problem to the problem of hitting a translate of an integral 
polyhedral cone represented by generators. To reduce this problem to the problem of 
hitting a rational cone we use some additional technical tricks and their description is 
contained in Subsection l5.3l As an intermediate step we consider the case of a simpli- 
cial cone (the generating vectors are linear independent). The general case is reduced to 
the simplicial one using representation of an arbitrary integral cone as a union of a finite 
set and a finite family of translates of simplicial integral cones (Theorem[6]). This result 
may be regarded as a sort of an integral Caratheodory's theorem. Recall that integral 
variants of the Caratheodory's theorem are actively studied (see, for instance, (4] |3) 
and have important applications to combinatorial optimization. 

Finally, in Section[6]we present some decidable and undecidable problems closely 
related to PB-realizability problem thus demonstrating its 'borderline' status with re- 
spect to computability. 

1 Algebraic and combinatorial properties of LRS 

Below we remind some algebraic and combinatorial properties of LRS which will be 
used in the sequel. All needed proofs and references may be found in ifMl l2l 171 . 

1 . Let A and h be, respectively, a linear operator and a linear function on a vector 
space V and let x, y G V. Then the generating function 

f A , x , y {t)=Y,KA k x-y)t r (2) 

is rational. 

2. The generating function of any LRS is rational and Taylor's series expansion 
coefficients of an arbitrary rational function at any point in its domain form a LRS. 

3. The set of LRS is closed under componentwise addition and multiplication (i.e. 
under Hadamard product of sequences). 

In the standard setting, the input to the Skolem problem is an integer LRS. But if we 
are interested in signs of the elements of LRS involved only, then there is no difference 
between integral and rational cases. 
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Indeed, let a LRS be specified by rational data and let N be the LCM of the de- 
nominators of all dj and 6 4 . Form an integral LRS 



Vn = J2 Ni aiVn-i ifn>d, y n = N n+1 b n if 1 < n < d. (3) 

i=l 

Proposition 1. N n+1 x n — y n foralln. 

Proof. For 1 < n < d the statement follows from the definition (fj). For n > d we 
proceed by induction. Assuming that the statement holds for all n < k, we get 

d d d 

y k = j2 Nia iVk-i = Y. N%a ^ xk - Nk ~ l+1 ) = J2 Nk+la i Xk -i = Nk+lx k- 

2—1 2—1 2—1 

Thus the statement holds for n = k. □ 
Remark 1. Trivially, using the Euclidean algorithm, the identity 

xy 



LCM(x,y) = 



LCF(x,y) 



and by Proposition [T] we can compute y n from x n in time polynomial in the length of 
the input data. So for algorithmic problems, which are concerned with the signs of LRS 
elements only, integral and rational cases are equivalent w.r.t. polynomial reductions. 

Hereinafter we assume that in algorithmic problems an LRS is represented by the 
list of coefficients a^, bj written in binary. 

Some LRS are solutions of enumeration problems and for this reason are nonneg- 
ative. For instance, the family of the so called N-rational sequences coincide with the 
set of generating functions of regular languages (see, S3). In the sequel N-rational 
sequences are called regular sequences. 

To be more precise, any deterministic automaton A over the alphabet S defines 
LRS s n , where 

Sn(A) = #{w : word w is accepted by A and \ w\ = n}. (4) 

Exactly the sequences of the type © will be called regular. 

LRS are not regular as they may have negative elements, but nevertheless the fol- 
lowing statement holds. 

Theorem 1 ( Mill Cor. 8.2]). Any LRS is a difference of two regular sequences. 



2 Orbits of linear maps and LRS 

In this section we state formally the orbit description problem and the related chamber 
hitting problem, and indicate relations of these problems to LRS. 
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Let hi, ... , h m be a family of affine functions defined on a coordinate space Q d . 
The sign patterns of these functions induce a partition of Q d into chambers. Formally, 
let s € {±1, 0} m . Then a chamber is a set 

{x e Q d : sign(hi(x)) = s» for 1 < i < m}, 

where sign(i) is a standard sign function 

1, fori>0, 
sign(t) = { 0, for t = 0, 
-1. fori<0. 



Orbit description problem (ODP). 

INPUT: a square matrix $ of order d; e?-dimensional vector xo\ a family of affine 
functions hi, ... , h m on Q d and a set of sign patterns si, . . . ,s r € {±1, 0} m . 

OUTPUT: 'yes' if any point of the orbit Orb$ xq falls into the union of chambers 
H Sl U, . . . , UH Sr , and 'no' otherwise. 

Remark 2. We assume that matrices, vectors and affine functions in the ODP and all 
related problems are represented by the component lists, where the components are 
written in binary. 

Chamber hitting problem (CHP) has the same input as the ODP, but there is one 
sign pattern s only. 

The output of CHP is 'yes' if the orbit Orb$ x intersects the chamber H s and 'no' 
otherwise. 

Proposition 2. CHP is Turing-equivalent to ODP. 

Proof. To reduce the CHP to the ODP we take the complement in the set {±1. 0} m to 
the sign pattern of the chamber in the input of an instance of the CHP and fix all other 
input parameters. We obtain an instance of the ODP that outputs 'yes' iff the instance 
of the CHP outputs 'no'. 

The reduction in the opposite direction is proved analogously. For any chamber not 
included in the input list of an instance of the ODP we solve the corresponding CHP. 
All the answers are 'no' iff the instance of the ODP reports 'yes'. □ 

Remark 3. It is rather obvious that the CHP is recursevely enumerable (belongs to 
the class S 1 of the arithmetic hierachy) and the ODP is co-enumerable (belongs to the 
class II 1 ). The proof of Proposition|2]shows that these problems are complementary in 
a broad sense. 

Another pair of complementary problems is the chamber description problem (the 
ODP restricted to the case of one chamber) and the problem of hitting the complement 
to a chamber. 

We do not know reducibilities between the chamber description problem and the 
CHP. 
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Proposition 3. The nonnegativity problem is equivalent to the ODP restricted to one 
linear function and a nonstrict inequality ( the chambers Hq, H + \). 

The Skolem problem is equivalent to the CHP restricted to one linear function and 
equality ( the chamber Hq). 

Proof. As is well known the LRS (Q]) is related to a linear operator A in the d-dimen- 
sional coordinate space given in matrix notation by 



.4 
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(5) 



\0 ... 1 0/ 
It is easy to check by induction that 

A k (b d , ...,&i) T = (x k +d,x k +d-i, ■ ■ ■ ,Xk+i) T '■ 

Thus, the Skolem problem for the LRS CD is reduced to the CHP with $ = A, 
%o — (bd, ■ ■ ■ , bi) T , hi — x d and the chamber Hq. 

Analogously, the nonnegativity problem is reduced to the ODP with the same input 
as above and two chambers H , H + \. 

The reductions in the opposite direction use standard facts about LRS listed in sec- 
tionQ] Indeed, the generating function h(A n xo)t n is rational. Hence, the sequence 
h(A n xo) is an LRS and both reductions follow. □ 

We don't know whether the general ODP and CHP can be reduced to, respectively, 
the nonnegativity problem and the Skolem problem, but we can indicate particular 
cases when such reductions do exist. 

Union of subspaces hitting problem (SHP). 

INPUT: a square matrix $ of order d; a d-dimensional vector a; ; a family of linear 
functions hjk on Q d , 1 < j < m, 1 < k < r,-. 

OUTPUT: 'yes' if the orbit Orb$ xq intersects the union of affine subspaces 

m 

[J{x : hji(x) = h o2 {x) = ■■■ = h jrj (x) = 0} 
i=i 

and 'no' otherwise. 

Lemma 1. The SHP is reducible to the Skolem problem. 

Proof. It follows from properties of LRS listed in Section[T]that the sequences 

VnC?) k) = h jk ($ n x ), j = 0, . . . ,m, 
are LRS. But as LRS are closed under componentwise sum and product, the sequence 

m Sj 

<Pn = n^2<Pn{j,kf 

3=0 k=l 
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is also LRS. 

Note that taking the input of the SHP we can algorithmically compute the repre- 
sentation of (p n in the form (|T). And the reduction of the SHP to the Skolem problem 
follows as (p n — iff for some j it holds 

<Pn(j> !) = Vntf, 2) = • • ' = (Pn(j, Sj) = 0. 

□ 

Remark 4. Note that m-reducibility in the proof above may not be a polynomial one 
as the degree of the sequence ip may be exponential w.r.t. the initial degree d. But the 
answer to a general SHP is disjunction of answers to separate subspace hitting prob- 
lems. And for this particular case the reducibility described in the proof is polynomial. 
Hence, we conclude that the SHP can be solved in polynomial time with the Skolem 
problem oracle. In other words, the SHP is polynomial time Turing reducible to the 
Skolem problem. 

The famous Skolem-Mahler-Lech theorem lf2ll7l [T5ll asserts that the set of zeroes of 
a LRS consists of a finite set and a union of a finite number of arithmetical progressions. 
The proof of Lemma[T]implies that the same is true for the SHP. 

Corollary 1 (Skolem-Mahler-Lech theorem for SHP). The set of those kfor which 
& k Xo falls into the union of subspaces Vi,... , V m , is a union of a finite set and a a 
finite number of arithmetical progressions. 

Polyhedron localization problem (PLP) is a particular case of ODP when the 
family of chambers forms a closed convex polyhedron (i.e. a solution set of a system 
of nonstrict linear inequalities). 

Lemma 2. The PLP is reducible to the nonnegativity problem. 

Proof. W.l.o.g. we assume that the sign pattern of the interior of a polyhedron is 
(+1, ■ • ■ , +1) and sign patterns for the faces of the polyhedron are obtained by re- 
placing some l's by 0's. Indeed, all negative components in a sign pattern can be 
removed by a sign change. And all zero components in the sign pattern of the interior 
of a polyhedron can be removed after checking that the orbit falls into (affine) subspace 
resulting from the solution of a system of linear equations (corresponding to zero com- 
ponents of the sign pattern). The last problem is certainly decidable as the affine hull 
of the orbit coincides with the affine hull of the first d + 1 orbit points. 

Now rearrange the sequences <f n {j) = hj($ n Xo) into one sequence <p n by con- 
sequently putting the elements of <p n (j) into places with an index having residue j 

def 

modulo m, i.e. <p S m+j — <Ps(j)i s > 0, < j < m. The resulting sequence <p n is an 
LRS. Indeed, its generating function f(t) is rational as it satisfies the identity 

m 

/(*) = E*w ro )' 

where for any j = 1 , . . . , m, a function fj (t) is a generating function for ip n (j ) . Hence, 
f(t) is a finite sum of rational functions. □ 
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The Skolem problem is the weakest of all problems mentioned above. To decide 
the Skolem problem it is enough to enumerate the family of LRS with nonnegative 
elements or the family of the PLP instances with 'yes' answers (due to m-reducibility 
in Lemma[2]). 

The following proposition belongs to folklore. 

Proposition 4. If the family of LRS with nonnegative coefficients is recursively enu- 
merable then the Skolem problem is decidable. 

Proof. If a sequence x n is a LRS then the sequence x\ — 1 is also a LRS as the family 
of all LRS is closed under componentwise sum/difference and product, see, SectionQ] 
But the sequence x 2 n — 1 has nonnegative elements iff x n ^ for all n. To complete 
the proof we apply E. Post theorem lfl3l . We enumerate all LRS with nonnegative 
coefficients and compare them to x\ — 1. In parallel we enumerate all LRS having 
zeroes and compare them to x 2 t — 1 also. One of these enumerations should eventually 
stop. □ 



3 CHP and regular languages 

In this section we relate the Skolem problem to some properties of regular languages. 
To be more exact we define problems of regular realizability: whether a given regular 
language contains a word of a particular kind. It is convenient to describe the problem 
of realizability as follows. Let L be a language in a finite alphabet E. Informally, 
L encodes a property that is checked. The input to the problem of L-realizability 
is a description of some regular language R C S* and we are asked whether the 
intersection L n R is nonempty. 

Remark 5. Depending on a way of presentation of a regular language we obtain dif- 
ferent in general variants of L-realizability problem. We assume that R is given via 
deterministic automaton accepting it. These clarification is of course immaterial if we 
are interested in decidability only, but may be important for complexity estimates. 

The permutation filter Pa is a language over the alphabet 0, 1} of all words 
#wi#w 2 # . . . wn#, where N = 2™, n > 1, and the set {w,}, i = 1,2, ... ,N 
consists of all binary words of the length n. 

Words from Pb are called permutation words. Informally, a permutaton word cor- 
responds to a permutation of a set of all binary words of some fixed length n. This 
length is called the block rank of a permutation word. 

Theorem 2. CHP and Pa-realizability problem are Turing equivalent. 

To prove the theorem we construct reductions in both directions, but as a matter of 
fact, their complexities are essentially different. 

In the next section we polynomially reduce the CHP to the Pb -realizability prob- 
lem and obtain as a corollary TVP-hardness of the last problem. This follows from 
the fact that the Skolem problem is iVP-hard [3] and according to proposition [3] it is 
polynomially reducible to CHP. 
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Remark 6. A different proof of TVP-hardness of the Pfl-realizability problem can be 
derived from the properties of periodic filters, see [ 18 1. 

The reduction in the opposite direction is shown in Section|5]but it takes superex- 
ponential time. 

4 Polynomial reduction of the CHP to 
the Pa-realizability problem 

We briefly sketch the idea of the reduction. First, using © we associate a LRS with 
any linear constraint (equality or inequality) in the description of a chamber. Then 
by proposition Q] we convert this LRS into an integral one. According to Theorem Q] 
this LRS can be represented as a difference of two regular sequences. It turns out 
that this construction is not enough to establish polynomial reduction and we elabo- 
rate it slightly and represent the integral LRS involved as a difference of two regular 
sequences on some explicit arithmetic progression of its indices. The last represen- 
tation can be computed in polynomial time. Next we find a special automaton that 
compares numbers of words of specified length in two regular languages provided the 
input is a permutation word. This trick gives a reduction of the Skolem problem to 
the Pn-realizability problem. To finish the proof we need an additional construction 
converting several automata described above into one. 
We proceed to the proof. 

At first note that for any square matrix A of order d, d-dimensional vector xo and 
a linear function h having rational coefficients there is a polynomial time algorithm to 
recover the coefficients of LRS for the sequence h(A n xo)- Indeed, the degree of the 
LRS involved does not exceed the order of the matrix and, hence, the coefficients of 
LRS and the initial data can be determined if we solve the corresponding system of 
linear equations for the first 2d elements of the sequence. 

In this section we assume that all LRS involved have integer coefficients and inte- 
gral initial data. This is possible due to Proposition [TJ as we are interested in signs of 
the expressions h(A n Xo) only. 

Next we express a LRS as a sum of weights of the walks in a digraph with fixed 
start and finish vertices. Let T(V, E) be a digraph (loops and parallel edges allowed) 
and let c : E —> Z be a weight function for edges of T. A weight of a walk is a product 
of weights over the edges of the walk. 

Lemma 3. There exists a polynomial time algorithm that takes an integer LRS {x n } 
and outputs a weighted digraph G x with two marked vertices s and f, such that the 
sum of weights over all walks of the length n starting at the vertex s and finishing at 
the vertex f equals x n for all n > 1. 

Proof. Take a LRS ([T]i and construct a digraph G x as follows. (See Fig.Q]) 

Connect the start vertex s and the finish vertex / by d vertex-disjoint (except for the 
terminals) paths Pi, . . . , P c i such that the edge length of Pi is i. Add d vertex-disjoint 
(except for the terminal /) directed cycles C\, . . . , Cd, passing through the vertex /. 
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Figure 1 : A graph for a LRS 



The edge length of C\ is also i. The cycles d have no common vertices with the paths 
Pj except for the vertex /. 

Now set edge weighting of G x . All edges of G x except for the first edges of all 
paths Pi and cycles Cj has weight 1. Set the weight of the first (outgoing from s) edge 
on the path Pj to Pi, i = 1, . . . , d, and set the weight of the first (outgoing from /) 
edge on the cycle Cj to Qi, i = 1, . . . , d. 

Let z n be the sum of weights over all walks of the length n from s to /. There is a 
recurrence on z n . Note that for n < d there exists a walk of the length n from s to / 
along a path of the digraph and for n > d any walk of the length n is obtained from a 
shorter one by adding a cycle. Thus 

i-1 

z i = Pi + ^ 1i-j z ji for 1 < z < d; 

d 

z n = ^2 <ljZ n -j, for n > d. 

To ensure equalities z n — x n set 

Qj = a 3- 

For 1 < i < d equalities Zi = b t imply 

Pi = h; 

P2 = h- qih; 

P3 =h- qxb 2 - q2h] 

To complete the proof note that all computations above take a polynomial time. □ 

Lemma 4. There exists a polynomial time algorithm that takes positive integers n ( in 
binary) and k (in unary) and outputs a digraph, a start vertex s and a finish vertex f 
such that the number of paths of the length k from s to f equals n provided k > log 2 n. 

Proof. Let L t be a doubled directed path having t+1 vertices {1, 2, . . . , t+1} and pairs 
of parallel edges (i, i + 1) between any consecutive vertices i and i + i = 1, . . . ,t. 
Clearly, there are exactly 2* walks of the length t between vertices 1 and 
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A thread L'j is a digraph obtained from Lj by attaching a directed path of the 
length k — j to its finish vertex j + 1. The start and the finish vertices of a thread are, 
respectively, the start vertex of Lj and the finish vertex of a path. 

Let {ji, . . . ,ji} be the positions of nonzero bits in the binary representation n. 
Take the threads Lj t , i = 1, . . . ,1 and identify their start vertices and, respectively, 
their finish vertices into global start vertex s and global finish vertex /. By construction 
there are exactly n paths of the length k from s to / in the resulting graph. 

The case of n = 11, k = 4 is illustrated in Fig. [2] 




Figure 2: Construction of LemmaHfor n = 11, k = 4 

Evidently the algorithm is polynomial. □ 

Now we are able to prove a variation of TheoremQ] 

Theorem 3. There exists a polynomial time algorithm that takes an integer LRS x n 
and outputs an integer i that is polynomially bounded in the input length ( i.e. the total 
binary length of all coefficients and the initial data of the LRS) and two deterministic 
automata A, B over the alphabet {0, 1} such that x n — S( n {A) — se n (B) for all n > 1. 

Proof. The first step is described in the proof of Lemma [3] As a result we obtain a 
digraph H . 

Then we choose integers M > log 2 max{|ai|, . . . , |a<j|, \bi\, . . . , \bd\} and k > 
log 2 (2Md). Now the integer I to output is set to M + k. 

At the next step we apply Lemma [4] and transform the digraph H to a digraph 
satisfying two additiional properties: 

1) the fan-out of each vertex is 2; 

2) the absolute value of each weight is 1. 

The transformation consists of the following steps. 

1. For each integer a; (6;) apply Lemma|4]to construct a digraph Ai (Bi) having 
exactly |et;| (\bi\) paths of the length M from the start to the finish. 

2. Replace the first edges of all cycles Ci (respectively, of all paths Pi) of the 
graph H by pasting in the graphs Ai (respectively, Bi). Namely, the start vertex of a 
graph Ai (Bi) is identified with the initial vertex of an edge and the finish vertex is 
identified with the end vertex of the edge. 

Other edges of the graph H are replaced by directed paths of the length I. 

Fan-outs of all vertices in the resulting digraph H' are at most 2 except for the start 
vertex s and for the finish vertex /. Note that fan-outs of the start and the finish vertices 
are at most 2Md. 

3. Attach to the start and to the finish vertex of the graph H' rooted directed binary 
trees of depth k having 2M d leaves. Replace the starting vertex of each edge outgoing 
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threads 



threads 




Figure 3: A graph H' 



the start vertex (the finish vertex) by a leaf of the corresponding tree in such a way that 
the fan-out of each vertex in the resulting graph H" is at most 2. 

As a result we obtain a modified digraph H" (see Fig. [3]). It has vertices with the 
fan-out less than 2. Add an auxiliary vertex u and edges to the u from all vertices with 
the fan-out less than 2 in the digraph H' to make all fan-outs equal 2. Then add two 
loops at the vertex u. After these operations we obtain a digraph G. 

Now we assign weights to the edges of the digraph G. The weight of an edge 
ingoing to the start vertex of a thread of a digraph A$ (Sj) coincides with the sign of 
the corresponding integer a* (6j). All other edges have weights 1. 

It is clear from the construction that in the digraph G the sum of weights over all 
walks of the length In from s to / equals x n . 

Now we construct automata A and B over the alphabet {0,1}. For both of them the 
state set is {±1} x V(G). Choose a labelling of the edges of the digraph G by and 
1 such that for each vertex edge labels of the two outgoing edges are different. (Recall 
that the fan-out of any vertex of G is 2.) 

Let's describe the automaton A. 

The initial state of the A is (+1, s) and the only accepting state is (+1, /). Reading 
a symbol a in the state (a, to) the automaton A goes to the state (a', to') if the edge 
(to, to') in the digraph G has label a. If the edge (to, to') has positive weight then 
a' = a. Otherwise, a' = —a. The rule is illustrated in Fig. [4] 

The automaton B differs from the automaton A in the accepting state only. The 
accepting state of B is (— 1, /). 

It is clear that the lengths of all words accepted by the automata A and B are 
multiples of I, 

Note that both automata can be constructed from the initial LRS in polynomial 
time. 

To finish the proof we note that x n equals the difference between the number of 
words of the length in accepted by the automata A and B. Indeed, the walks of negative 
weight correspond to the words accepted by the B while the walks of positive weight 
correspond to the words accepted by the A. □ 
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G {-1} x V{G) {+1} x V(G) 




Figure 4: Counting the sign 



Now we are ready to present a polynomial reduction of the Skolem problem to the 
i-fe-realizability problem. 

Theorem [5] implies that for the reduction it is sufficient to build up an algorithm 
running in polynomial time such that takes a pair of automata A, B over the alphabet 
{0, 1} and an integer £ in unary and outputs the description of an automaton C such 
that L(C) r\P B j= iff st n (A) = s (n {B) holds for some n. 

Let informally explain how to construct such an automaton. The automaton C 
is made of two automata C and C" by the product construction. The automaton C 
accepts iff C and C" accept. 

Recall that we are interested in operation of automata on permutation words, i.e., 
words of the type #w\#W2# ■ ■ ■ Wk#, where Wi £ {0, 1}*, i = 1, . . . , k. 

The block rank (the length of each wi) should be a multiple of I. This last condition 
is verified by the automaton C", It counts the length of wi modulo I and rejects in the 
case of a nonzero residue. Otherwise it accepts. 

The automaton C starts from an accepting state qo and reads the current block 
word Wi separated by the delimiters #. 

If Wi belongs to Lab = L(A) \ L(B) then C reads the next block word Wi+\ 
and checks whether it belongs to Lb a = L(B) \ L(A). If the last check fails then C 
rejects, otherwise it returns to the state go- 
lf Wi belongs to Lb a then C rejects. 

If Wi belongs to L^ = (L(A) ("1 L(B)) U L(A) U L(B) then C passes to a new 
accepting state q\. 

Starting from the state qi the automaton C reads the current block word Wi. If the 
block belongs Lab U Lb a then C rejects. Otherwise it returns to the state q\. 
The structure of the automaton C is pictured in Fig. [5] 

Let prove that C gives the required reduction of the Skolem problem to the Pb- 
realizability problem. 

First we must show that if an instance of the Skolem problem has a positive answer 
then C accepts at least one word from the permutation filter. 

Indeed, assume that some x n vanishes. Then si n (A) — si n (B) and there is a 
one-to-one correspondence between the words of the length In in Lab and the words 
of the length £n in Lb a- A required permutation word W is obtained by arranging 
corresponding words in pairs (other words of the length In are arranged arbitrary after 
all pairs). By construction C accepts W. 
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On the other hand, assume that C accepts some word W from the permutation 
filter. It follows from the construction that the block rank of W is a multiple of I. We 
can write W in the form W1W2. Here W\ is prefix of W that is a concatenation of 
#Wi#Wi+i, where Wi <E Lab and itfj+i £ Lba- And W2 is a suffix of W that is 
concatenation of where Wi € L^. Thus, |L(j4) \ = \L(B) \ L(A)\. It 

means that x n = and the instance of the Skolem problem has a positive answer. 

Clearly the size of the automata C is polynomial in the size of the LRS. This fin- 
ishes the proof of the reduction. 

From the results of this section and the results from [3 1 we get an immediate corol- 
lary. 

Theorem 4. The problem of P^-realizability is NP-hard. 

Now we extend the previous result to the general chamber hitting problem. 

The construction of the automaton C can be easily modified to accept a permutation 
word satisfying the condition sg n (A) < si n (B). In the notation above to be accepted 
the suffix W2 of a permutation word should contain at least one more occurence of a 
word from L(B) \ L{A). 

The counter part of the modified automaton C< does not change. The structure of 
the automaton C' K is pictured in Fig. [6] 

It is evident that the size of the modified automaton C< is upperbounded by the 
size of the automaton C up to a constant factor. 
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Figure 6: The structure of the automaton C<. The initial state is q . The only accepting 
state is q± 

To complete a reduction of the chamber hitting problem to the Pfl-realizability 
problem we construct an automaton that checks several conditions of the form 
Si n (A) = s in (B) (s en (A) < St n (B)). 

Repeating the previous argument let's construct for each constraint hi($ n xo) < 
(or hi(& n xo) — 0) from the description of the chamber a pair of automata Ai, Bi 
and an integer I common to all pairs Ai, Bi such that a constraint hi(Q n xo) < 
(hi(Q n xo) = 0) holds iff an inequality si n (A) < S£ n (B) (an equality sg n {A) = 
Shi(B)) holds. The main difficulty here stems from the condition that the block ranks 
in of all permutation words certifying the corresponding constraints in the reducibility 
must be equal. To solve this problem we use the following arguments. 

At first, choose the same integer i for all pairs of automata. The required value of I 
is 1 plus the three times the maximum of all logarithms of all data for all LRS involved. 
Then increase to £/2 the values of the parameters k and M from the proof TheoremfJ] 
It can be done by attaching paths to threads of the graphs constructed in Lemma|4] 

Next step is to construct for each pair Ai, Bi, where 1 < i < m and m is the 
number of constraints of CHP, the automaton C[ that certifies the inequality si n {A) < 
stn(B) (or the equality sg n (A) = sg n (B)) as described above. 

If m is not a power of 2 then by definition an automaton C[, where m < i < 2 P , 
p = [log 2 to], is a dummy automaton accepting all words. 

Now all automata C[ are combined into an automaton C that checks all conditions 
sin{A) ? si n (B) for some length in. 
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inCj: q' 



z represents neither i nor i — 1 

to the reject state to the initial state of Ki+\ 
Figure 7: Modification of a ^-transition in an automaton C[ 

The combined automaton has a product form C = C' x C" and accepts words 
#u>i#u>2# ■ ■ ■ wjv# from the permutation filter that satisfy the following require- 
ments. 

- The block rank is p + In for some n (in the sequel we denote Wi — UiVi, where 
\u\ = p). 

- Prefixes ui form a nondecreasing sequence w.r.t. the lexicographic order. 

- For each 1 < r < m let z be the binary representation of r — 1 of the length 
p. Form a subword #u>i#Wi+iif . . . #Wi+j# that consists of all block words 
with the prefix z. Then the word #Vi#Vi+i# . . . ifvi + j# is accepted by the 
automaton C' r . 

The first requirement is verified by a counter C" . The second and the third require- 
ments are verified by an automaton C , which is a sort of a concatenation of modified 
automata C[. In a modified automaton C[ each ^-transition from a state q' to a state 
q" is changed by a ^-transition from q' to the initial state of a copy of an auxiliary 
automaton Ki. This automaton reads a word z of the length p and compares it with 
p-length binary representations of i — 1 and i. If z represents i — 1 then Ki passes to 
the state q" of the automaton C[. If z represents i and q" is an accepting state in the 
automaton C[ then Ki passes to the initial state of the automaton C' i+1 . Otherwise it 
rejects. This modification is illustrated in Fig. [7] 

The initial state of the automaton C' is the initial state of C[ . The accepting states 
are the initial states of those copies of the auxiliary automaton that are inserted 
instead of ^-transitions to accepting states of the automaton C' 2P ■ 

The size of the automaton C is 0(Dm\ogm), where D is the maximum of the 
sizes over all automata Cj. Indeed, an auxiliary automaton Ki can be implemented 
using 0(p) = O(logm) states. So the size of a modified C[ is increased by a factor 
O(logm) and we combine sequentially 2 P = 0(m) modified automata C[. Let prove 
that the resulting automaton C gives the required reduction of the CHP to the Pb- 
realizability problem. 

Suppose that for some n all constraints in an instance of the CHP are satisfied. 
Then there exists a permutation of block words of the length In, which satisfies the 



Ki reads p-bit word z 



q" is rejecting 



z represents i 



q is accepting 
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accepting requirements of the automaton d for each 1 < i < m. Extending blocks by 
prefixes Zi that represent i — I'm binary and arranging permutations in the order of i 
we obtain the permutation word satisfying all three requirements stated above. Hence 
this word is accepted by C. 

On the other hand, assume that C accepts some word W from the permutation filter. 
The counter C" accepts. So the block rank of W is p + in. The automaton C accepts 
also. It follows from the construction that it must pass through all automata C[. There 
are 2 P possible prefixes of the length p. So the construction of modified automata 
implies that the sequence of block prefixes of the length p is uniquely determined. 
Namely, it has a form 

p bits p bits 

ooocl — o, ...,0000...0,... . 

V v ' 

2 ln prefixes 2 ln prefixes 

Thus acceptance by C' means that suffixes corresponding to each prefix z, where z is a 
binary representation of i — 1, form a permutation of the binary words of the length in 
and the corresponding permutation word is accepted by the automaton Cj. The latter 
implies that zth constraint in the instance of the CHP is satisfied. It implies that the 
answer for the instance of the CHP is positive. 




5 Reduction of the Pe-realizability problem 
to the CHP 

The reduction is composed from several intermediate reductions. 

At first, we show that the Pa-realizability problem is reducible to the Walk Weight 
Hitting Problem. The walks in this problem are walks in a digraph. 

Let r(V, E) be a digraph. We assume that the edges of the digraph are colored in s 
colors from the set {1,2,..., s}. 

For a walk r we define a weight w(t) 6 Z s as an s-dimensional integer vector 

W ( T ) = (ci(t),c 2 (t),...,c 1 {t),...c s (t)), 

where Ci (r) equals the number of edges colored in the color i on the walk r. 
Walk Weight Hitting Problem (WWHP) 

INPUT: a digraph T, an s-coloring of the edges of the T, two vertices a, b of T, an 
integer matrix $ of order s x s and an integer vector xq of dimension s. 

OUTPUT: 'yes' if the orbit Orb$ x intersects a set of weights of the walks from 
the vertex a to the vertex b and 'no' otherwise. 

Lemma 5. The P^-realizability problem is Turing reducible to the WWHP problem. 
The proof of Lemma|5]is presented in Subsection l5.ll 
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Then we will reduce the WWHP to the Integer cone Hitting Problem. An integer 
cone N(vi, . . . , v r ) is the set of vectors 

r 

2J CbiVi, a,; G N, 
i=i 

where Vi G Z d . We denote the set of nonnegative integers by N. 
Integer cone Hitting Problem (IHP) 

INPUT: a square matrix $ of order d; a d-dimensional vector xo\ a family of 
vectors vi G Z d , i = 0, 1, . . . , r. 

OUTPUT: 'yes' if the orbit Orb$ xq intersects the translate of the integer cone 
vq + N(vi, . . . , v r ) and 'no' otherwise. 

Lemma 6. The WWHP is Turing reducible to the IHP. 
The proof of Lemma|6]is presented in Subsection |5.2| 

Next step is to reduce the IHP to the Polyhedral cone Hitting Problem. This last 
intermediate problem is stated as follows. 
Polyhedral cone Hitting Problem (PHP) 

INPUT: a square matrix $ of order d; a d-dimensional vector xq; a family of linear 
functions hj on Q d ; a shift vector vq £ Q d . 

OUTPUT: 'yes' if the orbit Orb$ xq intersects the translate vq + K of the closed 
polyhedral cone 

K = {x G Q d : hj(x) > 0} 

and 'no' otherwise. 

The PHP problem is Turing-reducible to the ODP problem. The proof is similar to 
the proof of Proposition^ 

So, the total chain of reductions is 

Pu-realizability < T WWHP < T IHP < T PHP < T ODP < T CHP. 

As mentioned above the last two reductions are based on Proposition [2] So, to 
complete the reduction we should prove the following theorem. 

Theorem 5. The IHP problem is Turing reducible to the PHP problem. 

The proof of Theorem|5]is contained in Subsection |5.3l 

5.1 Proof of Lemma |5] 

We start from some preliminary work. 

Let i?be a regular language over the alphabet {0, 1, #} and let Abe a deterministic 
automaton with the state set Q accepting the language R. Our goal is to check whether 
R fl Pb 7^ 0. Let pass to the transition monoid and express this condition in terms of 
three maps /o, /i, /# of the set Q into itself induced by reading respective symbols. 
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Define f(w), where w — W1W2 ■ ■ ■ Wi £ {0, 1}, as 



lw e _ i+1 ■ 

1=1 



We denote the initial state of the automaton by q s and the accepting set by Q a . The 
reduction algorithm checks for each accepting state qf G Q a whether it can be reached 
from the initial state q s after reading some word from the permutation filter. Formally 
this condition means that for some word #w 1 #w 2 # • • • Wn# G Pa we have 

Qf = f#f( w N-*)^j ( 6 ) 

It is important that to verify the condition (|6]l we do not need to know the sequence 
Wi. It is sufficient to compute for each map g G the number v n (g) of its rep- 
resentations in the form f(w), where w <E {0, l} n (recall that N = 2"). Then the 
condition (O is rewritten as 



Qf 



= (jlf#9^f*q,, (7) 



where each map g G occurs exactly v n {g) times in the sequence gi. 

It turns out that the integers v n (g) can be expressed as the coordinates of the vectors 
taken from the orbit of a linear map. 

In Q-vector space Q(Q^) equipped with the basis {e(/)} indexed by maps f:Q—> 
Q we define a linear map by the action on the basis vectors 

< &e(/) = e(//o) + e(// 1 ). (8) 

(Recall that map composition is taken from the right to the left.) 

Proposition 5. v n {g) equals the e{g) -coordinate of the vector $™e(id) w.r.t. the basis 
{e(/ )}. Here id is the identity map. 

Proof. By induction on n. The case of n = is trivial. If the statement holds for 
n = k — 1 > then it holds for n = k: 

$ fc e(id) = $ Yl ^-i(g)e(g) = E $e(/M) = 

g£QQ ii>e{0,l} fc ~ 1 

= E (e(/N/o) + e(/H/i)= E (e(/(™0))+e(/M)) = 
we{o,i} k ~ 1 ive{os} k - 1 

= E </H) = E M9)e(g)- 

i«e{o,i} fc seQ Q 

□ 
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Now we are going to describe the condition (0 in terms of walk weights for a 
suitable graph. 

For this purpose we will use the Cayley graphs for monoids. 

Let G = {gi, . . . , g m } C be a set of maps. It generates a monoid M (a monoid 
operation is a map composition, recall that by definition a monoid contains also the 
identity map.). By definition the vertices of the Cayley graph Yq of the monoid are 
elements of the monoid M and (directed) edges have the form {h^gji) for h G M, 
gi G G. Note that the edges of the Cayley graph are naturally colored by elements of G. 
The Cayley graph of a monoid may contain parallel edges as the equality gih = gjh 
for i ^ j is possible for a monoid. 

Let Mqi be a semigroup generated by the maps fa, f\ of the automaton A. This 
semigroup is finite and can be described as the least set of maps from the state set Q to 
itself that contains the maps /o, /i and is closed w.r.t. multiplication by /o, /i. So, the 
semigroup Mqi can be constructed efficientljQ. 

In a similar way we define a monoid M generated by maps /#/, / G Mq\. This 
monoid is also constructed efficiently. 

Denote by Tm the Cayley graph of the monoid M w.r.t. the set of generators 
{/#/, / G M ,.\. 

It follows from the construction that 0) holds iff there exist an integer n and a 
walk t on the digraph Tm from the vertex id to a vertex h such that h(f#(q s )) = qj 
and {w{r))f # g — v n (g) holds for all g. And we conclude that the condition 01 is 
equivalent to the condition 

w(t) = $ n x . (9) 

Remark 7. Note that the size of the monoid M can be less than and $ acts on 

Q(Q^). To fix a difference in vector dimensions we extend the coordinates of walk 
weights by zero values for e(g) such that g ^ M, 

Thus the reduction algorithm solves with help of a WWHP-oracle several instances 
of the WWHP indexed by the accepting states qf and mappings h such that 

h(f#(q s )) = q f . 

An instance has the following input data: a digraph is the Cayley graph Tm, colors 
are generators, the initial vertex is the identity map id, the terminal vertex is h, the 
matrix is defined by (O and the initial vector is xq — e(id). 

If the answer is positive for some instance from the described set then the answer 
in the instance of the ffe-realizability problem involved is also positive due to Proposi- 
tion [5] and condition ||9). 

Otherwise, it follows from the definition of the graph Tm and Proposition [5] that 
the answer in the instance of the PB- rea li za bility problem is negative. 

Lemma|5]is proved. 

5.2 Proof of Lemma H] 

We need the following auxiliary lemma. 

1 Hereinafter 'efficiency' means a mere existence of an algorithm. 
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Lemma 7. Let a, b be vertices of a digraph T colored in s colors and let W(a, b) be 
the set of weights of all walks from the vertex a to the vertex b. The set W(a, b) is a 
finite union of translates of integer cones in If . 

The list of the cones and the vectors of translation is constructed efficiently. 

Proof. Writing a sequence of colors along a walk gives a word in the s-letter alphabet. 
The collection of such words for all walks from the vertex a to the vertex b forms a 
regular language. 

The weight of a walk is just the Parikh image of the corresponding word. So the 
statement of the lemma follows from Parikh's theorem [ 10 1. □ 

Lemma [7] implies that to check the condition (O it is sufficient to check several 
conditions of the form 



where vq in an integer s-dimensional vector and W is an integer cone in Z s . It gives a 
reduction of the WWHP to the IHP 

5.3 Proof of Theorem g] 

Let Wq C Q s be a rational cone generated by vectors V\, . . . , v r . There is an algorithm 
(see, e.g. IT21 §7.2]) for passing to dual polyhedral description of Wq. 

The cone Wq may contain some additional integral points not belonging to integral 
cone W so that conditions ( fTOb should be modified. 

Let us start with a simple case of an integral simplicial cone, when generating 
vectors Vi are linear independent. 

Proposition 6. vectors v\, . .. , v r are linear independent then representation x = 
diVi, Oj 6 N holds if and only if the following two representations hold: x — 
J2i hVi, bi G Q+ and x = £V c^, e Z. 

Proof. The implication is trivial and the opposite implication => follows from linear 
independence of vectors Vi, as the equality 



Vector x can be represented in a form x = J^i c i v i, where q e Z, iff x is an 
element of the subgroup G generated by vectors vi in Z s . And this condition in turn 
is a conjunction of two conditions. The first condition is simple: a; is in a subspace 
generated by vectors Vi in the vector space space Q s . To formulate the second condition 
let extend the set Vi of generators of G to the basis of Z s by adding some unit coordinate 
vectors ej so that the resulting system formed from Vi and ej constitute a basis of Z s . 
Such extension is always possible in view of linear independence of w^. Let G be the 
group generated by Vi and ej. 



<$> n x G w + W, 



(10) 




holds iff bi 



c-it i — 1 , . . . , r. 



□ 
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Proposition 7. Vector x belongs to the subgroup G generated by vectors in the 
group I s iff x is in the subspace generated by vectors Vi in Q s and x belongs to the 
subgroup G. 

Proof. The implication <= is obvious. To prove the opposite implication => assume 
that x G G. By construction of the subgroup G it follows that 

j 

where g 6 G and ej denote those added coordinate unit vectors. Now it follows from 
linear independence of Vi and ej that if vector x belongs to the subspace generated by 
vectors Vi in the coordinate vector space Q s then all coordinates Xj are zero. Hence x 
belongs to the subgroup G. □ 

It follows from the propositions|6]and|2]that a vector falls into integral conic hull of 
vectors if and only if three conditions are fulfilled. The first two conditions are Q- 
linear and control whether a vector falls into the cone Wq. The third condition consists 
in checking that a vector belongs to a fulldimensional lattice in Z s (a subgroup of V 
having finite index in Z s ). Checking the first and the second condition for the orbit 
of a linear map are particular variants of CHP and the third condition holds for a nice 
family of orbit points. 

Proposition 8. Let G = (i>i, . . . ,v s ) C V be rank s subgroup of 'I s , let vq be an 
integral vector in 1/ and let <£> be an integral s x s matrix. 
Then the set 

H = {n: <S> n x G v + G} 

is a union of a finite set Hq and a finite set of arithmetic progressions. Moreover, there 
is an algorithm to compute H. 

Proof. First let G be a diagonal subgroup whose generators are columns of the matrix 







. 


• °\ 


D = 







. 






. 


• Is) 



The shift induced by G is given by 

Xi = £i (mod qi), l<i<s. (11) 

Let q be the least common multiple of g,. Now calculating Q k xo modulo q we 
can find a finite set of exponents k, and the corresponding arithmetic progressions k 
(mod q), satisfying (fTTT i. Besides these arithmetic projections H may include only 
numbers that should be less than initial terms of the progressions found on the previous 
step, so that we can directly compute this finite set. 

The general case is reduced to a diagonal one by reducing the generator matrix M 
to the Smith normal form (see [ 12 , §4.4]) If G is a subgroup of Z s of full rank then 
Smith normal form is given by 

M = UDV, (12) 
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where U, V are unimodular matrices and D is a nondegenerate diagonal matrix. V 
corresponds to column transformations and gives a system of new generators v\, . . . , 
v 8 of G. U corresponds to row transformations and shows how the vectors of the 
initial basis can be expressed as a linear combinations of the new basis vectors ec 

(ei,...,e s ) = fa,... e a )U. 

By (fT2l the generator matrix Vi of the group G is diagonal in the basis e^. 

To complete a reduction to the diagonal case we should express the matrix $ and 
the initial vectors xq, vq in the new basis gj. □ 

Lemma 8. If vectors v±, . . . , v r are linear independent then the IHP is reduced to the 
PHP. 

Proof. Let us describe an algorithm with PHP-oracle that solves the IHP for instances 
of linear independent set of vectors v\, . . . , v r . 

First let find the set H defined in proposition [8] We should check whether con- 
dition ( [Tol l holds for H. The finite part Ho C H can be checked directly. And for 
each arithmetic progression n = no + Nk, k = 0, 1, . . . in H we compute a vector 
x\ = $™° and a matrix <£>i = $> N . Then using PHP-oracle we check whether there 
exists k such that $\xi — vq hits the rational cone Wq. If all such tests fail we answer 
'no'. Otherwise we answer 'yes'. 

To check correctness of the algorithm involved note that it follows from the propo- 
sitions [7] and [8] that all tests fail if and only if the orbit $"xo either does not intersect 
the cone Wq at all or hits it for exponents n such that $™xo — «o is n °t in the group G. 
Hence, the answer in the IHP is negative. 

On the other hand, if some test is successful then we find n, such that $™a:o — Vq is 
in the group G and hits the cone Wq. It follows from proposition|6]that the answer in 
the IHP is positive. □ 

Now we discuss the general case. Let 

W N = N(v l7 ... : v r ) 

be the set of integral linear combinations of v% with nonnegative coefficients (the in- 
tegral cone generated by Vi). Below we describe an algorithm that given any integral 
cone Wfi produces its representation as a union of a finite set Wo of singular points and 
a finite set of translates of integral simplicial cones m + N(v^ , . . . , Vi t ). 

Recall that Caratheodory's theorem (see, e.g. Ifl2l §7.7]) states that any vector in 
Q+(vx, . . . ,v r ) is a rational nonnegative linear combination of some < s (where s is 
the dimension) vectors from the system Vi, ...,v r . 

As one can form not more than Pj systems of linear independent vectors out of r 
vectors, then it is sufficient to consider the case of the intersection of the integral cone 
Wjh and a rational simplicial cone with s generators chosen from the set Vi, . . . ,v r . 

Let K be one of such simplicial cones and let v\, . . . , Vt be the set of its generating 
vectors. K contains N(v\, . . . , vt), but besides that K may contain some extra integral 
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points. Now we use the fact that K has Hilbert basis: a set of integral vectors u\, ■ ■ ■ , 
u m , such that 

K C\Z S =N( Ul ,...,u m ). (13) 
In particular, it follows from ( [T3l that 

tn 

k nz s - |J (ui + nivx,...,^)). 

Recall that one can effectively find Hilbert basis for K as it coincides with the set 
of integral points of the polytope 

{\ x vi + X 2 v 2 + ■■■ + X t v t : < Xi < 1}. 

Some vectors from Hilbert basis belongs to the subgroup Z s , generated by Vi, i.e. 
to the integral hull of u,;. Checking this condition is effective as it is reduced to solving 
linear Diophantine equations. 

Let take one of such vectors u = J2i^i v i an d l et B be maximum of modules of the 
coefficients hi. Now all points of the intersection of u + N(5i, . . . ,vt) and Wn belong 
to the union of the finite set 

t 

{x : x = u + aw, < di < B}, (14) 

i=i 

the integral cone 

u+CB + lJ^Wi + Nfa,...,^) (15) 

i 

and (B + l)t sets of the form 

(u + aVi + Q + (vi,V2, ■ ■ .,Vi-i,v i+1 , . . . ,v t )) nWfsj, < a < B, 1 <i < t. (16) 

Each set in ( fT6] l is an intersection of the integral cone Wm and a translate of some 
rational simplicial cone of the lower dimension t — 1. 

Now we can give a complete description of the algorithm that given an integral cone 
Wn and a translate of some rational simplicial cone whose generating vectors are taken 
from the set of generators of Wn, finds representation of the intersection of the cones 
involved as a union of a finite set of singular points W$ and a finite set of translates of 
simplicial integral cones. 

The procedure is recursive. Applying it to the intersection of and a i-dimen- 
sional rational simplicial cone K we at first compute Hilbert basis of K. Then for any 
vector in Hilbert basis that belongs to the integral hull of vectors v\, . . . , v r , we find 
the sets (fl4l . (TT3T > and ( [T6j >. 

The set (fl4l is added to the singular part Wo. The simplicial integral cone ( fTBI l 
is included to the family of integral cones computed on the previous steps and the 
algorithm is continued recursively by processing all sets (fTBI l. 

The procedure is finite as the sets dT6b are empty for any one-dimensional cone. 
Hence, the recursion tree has finite height and finite degree at any point. 

Thus the following theorem is proved. 
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Theorem 6. An integral translate of an integral cone 

v +N(«i, ...,v r ) 

can be represented as a union of a finite set and a finite family of translates of simplicial 
integral cones. There is an algorithm that finds such a representation from the list of 
vectors Vq,V\,..., v r . 

Now the main result of this section easily follows. 

Proof of theorem^ Using algorithm from Theorem [6] find representation of the inte- 
ger cone as a union of a finite set of singular points and finite family of translates of 
simplicial integral cones. Now for any cone in the family check whether the orbit hits 
it using PHP-oracle from Lemma|8] 

Checking whether or not the orbit hits a point (and hence, any finite set of points) is 
also reduced to the PHP as a point may be regarded as a translate of the zero cone. □ 

6 Decidable and undecidable variants of the regular re- 
adability problem 

The Skolem problem is open for almost eighty years. Using slight abuse of language, 
presently it falls 'on the border between decidability and undecidability' [7 1. In this 
section we show that analogous 'borderline' pattern holds for a more general P\- 
realizability problem and give some decidable and undecidable problems closely re- 
lated to it. All these problems are problems of regular realizability. The languages 
specifying the problems consist of binary block words separated by the delimiter 
All blocks have the same length (the block rank n). Blocks of a block word form a 
multiset of binary words of the same length. We will define languages indicating the 
properties of this block multiset. 

Now we present decidable examples. 

The surjective filter 5b consists of those block words which block multiset contains 
all words of the length n. 

The injective filter Zjg consists of those block words which block multiset is a set, 
i.e. each block appears at most once in a word from 2b- 

It is clear from the definitions that fl 5b = Pa- 

The problems of Is-realizability and Si-realizability are decidable. Let's outline 
the proof. It turns out that both 7r- or, respectively, 5b- realizability can be reduced to 
some restricted versions of the integer cone hitting problem whose decidability follows 
from specific properties of maps $ defined in Subsection l5.ll 

We use the componentwise partial order on the integer orthant 

x <y ^ Xi<yifor all i. (17) 

Up-hitting Problem 

INPUT: a square matrix $ of order d; a e?-dimensional vector xo\ a family of 
vectors Vi £ Z d , i = 0, 1, . . . , r. 
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OUTPUT: 'yes' if the orbit up-shadow intersects the translate of the integer cone 
vq + N(vi, . . . , v r ) — vo + W and 'no' otherwise. It means that 

$™x < y (18) 

holds for some integer n and y G i>o + W. 

The Down-hitting problem is defined similarly except for the condition ( fT8l which 
is replaced by the condition 

y ~< ® n xo- (19) 



Lemma 9. The I^-realizability problem is Turing reducible to the down-hitting prob- 
lem. 

Proof. Repeat the arguments from the proof of Lemma|5] Now it is possible that some 
binary words in a block word are missed. It means that the condition (0 is replaced by 

w(t) -< $"x . (20) 

Applying the arguments from Lemma [6] we see that (l20b is transformed to ( fT9l for 
cones appeared in the reduction from Lemma[6] □ 

In a similar way we reduce the surjective filter. 

Lemma 10. The Sa-realizability problem is Turing reducible to the up-hitting problem. 

Proof. In a block word taken from the Sb all binary words of the length n appear (pos- 
sibly, several times). It leads to the condition ( fT8l for cones appeared in the reduction 
from Lemma[6] □ 

Dickson's lemma claims that there are no infinite antichains in the poset (Zl, -<). 
So an orbit up-shadow is a finite union of translated copies of the orthant Zl. In the 
case of an orbit down-shadow copies of the orthant are replaced by 'parallelepipedons' 
(the Cartesian products of segments). Thus the up-hitting problem as well as the down- 
hitting problem is reduced to the nonemptiness check for intersections of integer cones 
(parallelepipedons), which is an integer linear programming problem, provided the 
representations mentioned above can be constructed efficiently. 

To construct the aforementioned representations we use specific properties of asymp- 
totic behavior of the orbit points $> n xo, where $ is determined by dS). 

Recall that an e(g) -component of $ n xo is expressed as the number v n (g) of walks 
of the length n from the vertex id to the vertex g in the graph T. The vertex set of the 
graph r is V(T) — and the edge set E(T) consists of pairs in the form (/, //o) or 
(/, ffi). (See Subsection l5.1l for a detailed exposition.) 

Note that the integer v n (g) is the number of words of the length n in a regular 
language. This language is recognized by an automaton with the transition graph T. 
Thus the generating function 

oo 

¥>»(*) = $>»(»)*" 

n=0 
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for the language is a rational function and its representation in the form P(t) /Q(t) can 
be found efficiently. 

In the arguments below we need some properties specific to generating functions 
of regular languages. So it is more suitable for our purposes to analyze the asymptotic 
behavior in combinatorial settings by considering walks on the graph V. 

Proposition 9. For any vertex g £ V(T) the set 

P g = {n: u n (g) > 0} 

is a semilinear set (a finite union of an exceptional finite set and finite collection of 
arithmetic progressions) and its description can be constructed efficiently. 

Proof. Regard the Y as the transition graph of a nondeterministic automaton over a 
1 -letter alphabet. Then the proposition follows from Parikh's theorem. □ 

Remark 8. Differences of all progressions in Proposition [9] are the cycle lengths in 
the graph T. So they are divisors of the least common multiple of integers from 1 to 
|V(r)|. In the sequel use we denote this common multiple by N. 

Take a vertex g on a directed cycle of the length I. The following inequality 

v„+i{g) > v n (g). (21) 

holds. Indeed, one can extend any walk of the length n by the cycle. 
Now we divide the vertices of the graph Y into three groups. 

- Vi consists of vertices v such that some directed cycle (possibly, a loop) goes 
through the vertex v. 

- V 2 consists of vertices v such that there is a walk starting at the id, finishing at 
the v and passing through a vertex from the set V\ . 

- V3 consists of all other vertices. 

Proposition 10. If g e V 3 then v n (g) = 0forn> \V(Y)\. 

Proof. Any walk of the length n > \V(Y) | from id to g must contain repeating vertices. 
It means that a part of the walk is a directed cycle. So the walk passes a vertex from 
the set Vi. □ 

Proposition 11. The inequality 

v n +N(g) > v n (g), (22) 
where N = LCM(1, . . . , |V(r)|, holds for all g e V(Y) andn> \V(Y)\. 

Proof. For g £ V3 apply PropositionfTOl 

For g G V\ the inequality (1221 follows from d2TT > and Remark|8] 

Now take a vertex g E V 2 . The set V g C V\ consists of vertices g' 6 V\ such that 

g' belongs to a walk from id to g and for all walks of this type all vertices after the g' 
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along a walk are in the set V%. The subgraph T g is induced by the edges of all walks 
from a vertex in V g to g. 

Let observe the following properties of the T g . 

There are no edges to vertices of the set V g in the graph T g . Indeed, such an edge 
contradicts definition of the set V g . 

There are no edges outgoing from the vertex g in the graph T g . Otherwise one 
would detect a directed cycle passing through g. 

From these properties we conclude that the graph T g is acyclic. By definition there 
are no directed cycles passing through vertices in the set V2. Other vertices in the T g 
are in the set V g . There are no directed cycles passing through these vertices because 
there are no edges ingoing to them. 

Note also that the maximum of the path length from g' 6 V g to g does not exceed 
\V(T)\. 

From all the properties above we get 

v n(g) = ^2 Pg',k v n-k{g'), (23) 
9' ev g 
fc<|v(r)| 

where p g ',k is the number of paths from g' and g in the T g . 

Applying the inequality d22l to all terms in the right-hand side of (1231 we get the 
same inequality for the left-hand side, i.e. for the vertex g. □ 

Theorem 7. The S^-realizability problem is decidable. 

Proof. It follows from Lemma [10] that it is enough to construct an integer cone rep- 
resentation for an orbit up-shadow. Proposition Ql] implies that the orbit up-shadow 
is 

N 
i=0 

where to is a dimension, i.e. the cardinality of Q®. 

So the problem is reduced to the integer linear programming problem. □ 

To prove decidability of the injective filter we should determine for each < r < 
N unbounded components of § nN+r XQ and the limit values of the bounded compo- 
nents. 

All subsequences t r g (n) = is n N+r(g), where < r < N, are nondecreasing due to 
PropositionfTTI 

A subsequence t g (n), where g G V3 stabilizes for n > \V(T)\ and the limit value 
for it is 0. 

It follows from d23l that a subsequence t r g {n), where g S V2, tends to infinity iff at 
least one of the subsequences t r g 7 (n) tends to infinity, where p g >.k 7^ 0. 

The remaining case t g (n), where g € V\, is covered by the following proposition. 

Proposition 12. Let g s V\. Then lim t g (n) = 00 iff there exist a directed cycle C 
passing through g and an edge (g' , g") such that 
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(i) the edge (g 1 , g") is not included in the cycle C; 

(ii) the cycle C passes through the vertex g"; 

(iii) v n N+r-i-i{g') > 0, where t is the distance from g" to g along the cycle C. 

Note that due to Remark [8] the conditions (iii) are equivalent for all n. It is clear 
that (i)-(iii) are verified efficiently. 

Proof of Proposition^^ 'If part of the proposition follows from 

V{n+\)N+r{d) > VnN+r{g) + V(n+l)N+r-t-l{g') > V-nN+r{g)- (24) 

Prove now 'only if. Suppose that v n N+r(g) = T for n > uq. For any directed 
cycle C passing through g and any edge satisfying (i)-(ii) the first inequality in d24l i 
implies that v n N+r-i-i{g') = for n > no. □ 

Theorem 8. The I^-realizability problem is decidable. 

Proof. Determine all unbounded components for all subsequences & nN+r XQ and the 
limit values for bounded components. For this purpose use Proposition [12] and the 
observations preceding it. Then the down-shadow is the union of the sets 

y g >0, if lim ($ nfi+r x )g = oo, 

n—toc 

y?>y g >Q, if Mm{* nN + r x ) g = vf. 

over all < r < N. Here y g are coordinates in the space Q® Q . 

So the problem is reduced to the integer linear programming problem. □ 

Now we present an undecidable realizability problem that is related to the P«- 
realizability problem. 

In the construction we use a track product of languages consisting of block words 
(block languages). Here blocks are words over a finite alphabet £ and a block word 
consists of blocks of the same length separated by the delimiter 

Let Si, £2 be finite alphabets. For the alphabet Si x S2 there are two natural 
projections from (Si x S 2 )* to S* (resp. to S|): 

7Ti : (ai, 61) (02, 62) • • • (a n , K) i-> a-ia 2 ...a n , ^ 
7T2 : {ax, bi)(a 2 , b 2 ) ■ ■ ■ {a n , b n ) ^ bib 2 ■ ■ . b n . 

The track product of two block languages L\ (over an alphabet Si) and L 2 (over 
the alphabet S2) is a block language L = Li\\L 2 consisting of all block words over 
the alphabet {#} U Si x E 2 such that the projection wi is in the language L\ and 
the projection tt 2 is in the language L 2 . (For consistency we assume that ~K\{jf) = 

7T 2 (#)=#.) 

Denote by Perg the block language consisting of all periodic words over the alpha- 
bet S (all blocks of a word in Pcrs are equal) and by the block language consisting 
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of permutation block words over the alphabet £ (blocks of a word in Ps form the set 
of all words in S n , where n is the block rank). 

The Pers-realizability problem is decidable. Actually it is PSPACE-complete 1 18 1. 
It turns out that the track product of the periodic filter with the permutation one is 
undecidable. 

Theorem 9. There are alphabets Si, £2 such that the (Perj^ ||Ps 2 )-realizability prob- 
lem is undecidable. 

A suitable undecidable problem that is reduced to the (Perj^ ||Pe 2 )-realizability 
problem is the following. 

Zero in the Upper Right Corner Problem. (The ZURC problem.) For a given 
collection of D x D integer matrices Ai, . . . , An check whether the multiplicative 
semigroup generated by {A;} contains a matrix M such that M%d = 0. 

In other words the ZURC problem is to check an existence of an integer sequence 
ji, . . . , je, where 1 < jt < N, such that 

{A n A n ...A u ) lD = Q. (26) 

Theorem 10 (see JT]). The ZURC problem is undecidable for N — 2 and D = 18. 

We will reduce the ZURC problem with N = 2 and D = 18 to the (Pcr Sl ||Ps 2 )- 
realizability problem where 

S 1 = [l,...,iV], £ 2 =[1,...,£>]x[1,...,£>]x{0,1}. (27) 

The reduction is similar to the reduction in Section [4] 

Let Ax, . . . , Ajy be an instance of the ZURC problem for D = 18, N = 2. Rewrite 
the matrices in the form 



/ j 3 j 3 3 3 \ 

I £ ll m ll £ 12 m 12 • ' • £ lD m lD \ 

3 3 j 3 3 3 

£ 21 m 21 £ 22 m 21 ■ ■ ■ e 2D m 2D 

, j j j j j j 

\£ D1 m D1 ... £ D (D-l) m D (D~l) £ DD m DD/ 



where m 3 ik > and e\ k € {±1,0}. Let M be the maximum of m\ k . 

Fix now a sequence Aj 1 , Aj 2 , . . . , Aj e . Matrix elements in the product have the 
form 

(A n A n . ..A h ). k = £e(r)m(r), (28) 

T 

where r runs over all sequences of pairs (i a k a ) such that the length of a sequence is I 

and i\ — i, kf — k, i a +\ — k a , 

< T ) = n 4: ka . >»(r) = n m t^ ■ ( 29 > 

01=1 a — 1 



30 



Using the expansion (l28l [29b we define the partition of words of the length £ ■ 
|~log 2 M] over the alphabet Si x S 2 into three sets T^(j), T^(j) and T^ k A {j), where 
j = ji, ■ ■ ■ ,je- (Below we drop out j while the sequence j is fixed.) 

It is convenient to represent a word over the alphabet Si x £2, where are given 
by (|27l i, by a 4-row table. A table column represents a symbol in the word. The first 
row bears symbols from Si while the remaining three columns represent symbols from 
S2 (they have three components as indicated in (l27Ti). 

A word in the set T+ (T~) can be divided in £ subwords of the length |~log 2 M~\ . 
The ath subword has the form 



f ja 


j a ■ 


ja \ 




* a 


z a 






k a 


\/?0 


/3i • 


■ P\log 2 M]-lJ 



(30) 



in the table representation defined above. The upper three elements in each row are the 
same in (13 Oi l. Note that j a is the ath element of the sequence j. The fourth row is a 
binary representation of an integer (3. 

A word in the set TX should satisfy the following requirements 

(a) ii = i; 

(b) k e = k; 

(c) i a+ i = k a ; 

(d) ft < ml2kj 

(e) e(r) — 1, where r is the sequence of pairs (i a , k a ) and e(r) is given by d29l . 

A word in the set T^ k should satisfy the requirements (a)-(d) and the modified 
requirement (e') e(r) = —1. 

The set T^ d collects the rest of words. 

Proposition 13. In notation above 

(Aj x Aj 2 . . -A n ). k = \T ik \ - \T ik \. 

Proof. Restricting words in the set TVt to the upper three rows one obtains all correct 
sequences r such that e(r) = 1. The same restriction for the set T.T gives the sequences 
r such that e(r) = —1. 

The multiplicity of a sequence r in the set TV£ depends on fourth rows of the ta- 
bles. According to the requirement (d) and the permutation property there are exactly 
m^° fc subwords bearing (j a , i a , k a ) in the upper three rows. So the multiplicity of the 
sequence r is m(r), where m(r) is given by ( 129) . □ 

It follows from the above construction that the sets Tf k (j), where $ e {+, — ,bad}, 
do not intersect for different j because the sequence j is recovered from the first row 
in the table representation. 
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Proposition 14. For a fixed collection of matrices A\, . . . , An and for each pair i, k 
the language 

T£ S = {jT? k (j) 

3 

is regular for $ £ {+, —.bad}. 

Proof. The requirements (a)-(d) are verified by local check on the subwords of the 
fixed length [log 2 M~\ . 

Computing the sign e(t) can be done in the cyclic group of two elements. 

So the sets T^' ± are regular. But the class of regular languages is closed under the 
complement. Thus the set T ik ' hdd is also regular. □ 

Proof of Theorem^ We repeat the construction from Section [4] The reduction algo- 
rithm takes an instance of the ZURC problem (N — 2, D — 18) and constructs an 
automaton C such that the condition d26l i holds for some element in the semigroup 
generated by the input matrices if and only if L(C) D Per^j ||Pe 2 - 

The automaton expects an input w from the language Perj^ \ \Ps 2 sucn that block 
rank is I ■ [log 2 M~\ and the word w is a certificate for zero representation d26i >. A 
period in the first (periodic) component determines the sequence j = j\ , . . . , ji such 
that d26l > holds. The automaton expects that each symbol in the sequence j is repeated 
[log 2 M~\ times. In the second (permutation) component the automaton expects that 
the blocks from the sets T^ D (j), T^~ D (j) are paired and are followed by the blocks 
from the set T^. Note that such a pairing exists iff \T^ D \ — \T^ D \. The structure of 
this part of the automaton C is similar to the automaton C" described in Section [4] and 
shown in Fig. [5] 

The correctness of the reduction is proved in a way similar to the arguments in 
Section[4] If the automaton accepts a word w in the language Per^ \\Ps 2 then one can 
extract the sequence j from the first components of symbols in the word w. The check 
in the second components guarantees that d26b holds for the sequence j. We apply here 
Proposition[T3l 

In other direction, if d26l i holds for a sequence j then there exists a word in the lan- 
guage Per^ ||Ps 2 satisfying the properties expected (and verified) by the automaton C. 
Thus L{C) n Per El ||Jfc 2 ^ 0. ' ' □ 

Remark 9. By suitable encoding of the symbols of the alphabets Si and S 2 one can 
prove that the (Peri||PB)-realizability problem is also undecidable. 
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